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The KadanofF-Baym (KB) equations are solved numerically for infinite nuclear matter. In particular 
we calculate correlation energies and correlation times. Approximating the Green's functions in the 
KB collision kernel by the free Green's functions the Levinson equation is obtained. This approxi- 
mation is valid for weak interactions and/or low densities. It relates to the extended quasi-classical 
approximation for the spectral function. Comparing the Levinson, Born and KB calculations allows 
for an estimate of higher order spectral corrections to the correlations. A decrease in binding energy 
is reported due to spectral correlations and off-shell parts in the reduced density matrix. 



I. INTRODUCTION 

The quantum Kadanoff-Baym equations (KB) Q de- 
scribe the time-evolution of the two-time (one-particle) 
Green's functions G{p,t,t'). Imposing various approxi- 
mations they have played an important role in the past 
developing corrections to the classical Boltzmann equa- 
tion such as memory-effect and damping. With some 
restrictions it is now however feasible to solve these equa- 
tions numerically without approximations. 

We like to emphasize the following differences between 
the Boltzmann and KB-equations. In the quasi - clas- 
sical Markovian Boltzmann-equation correlations are ex- 
plicitly neglected. The spectral functions are given by 
the quasi-particle approximation. The kinetic energy is 
conserved in each binary collision, while the correlation 
energy is neglected. The equilibrium density-distribution 
(in momentum space) is a Fermi-distribution. Oppo- 
sitely, in the KB-equations correlations are carried by 
the two-time Green's functions. The total energy includ- 
ing the interaction potential energy is conserved. The re- 
duced density distribution is given by a frequency integral 
of a spectral-function of non-zero width and a distribu- 
tion function. At equilibrium this distribution-function 
is again given by a Fermi-distribution. 

The various approximations of the KB-equations differ 
essentially by the reduction schemes of the two - time 
Green's function to the reduced density matrix or even 
to the quasiparticle distribution. That this reduction is 
nontrivial can be traced to the fact that the Fermi and 
quasi-particle distributions decrease exponentially with 
energy while the reduced density matrix possesses power 
tails. In this paper we address some common schemes 
from the numerical point of view. 

Numerical results of the quantum KB-equations al- 
ready have been compared in the past with the classical 
Markovian dynamics as well as with other frequently used 
approximations ||,^. The study was made in particular 



with reference to collisions between heavy ions usually 
studied with the Boltzmann-equation or the BUU,VUU 
etc versions thereof. The general conclusion of the study 
was that quantum-effects reduce relaxation-rates by a 
factor as small as one half. This fact can be interpreted 
as being due to the pole renormalization or wave function 
renormalization 

Since the first numerical applications of the KB- 
equations by Danielewicz several contributions to this 
evolving new field have been published with applications 
to nuclear matter to one- and two-band semi- 

conducto rs [^ ,p|,to phonon-production in e-e collisions in 
pla smas |10|| as well as to electron plasmas in general 
[|lT| , ^ . Details of the computational methods is pub- 
lished in Computer Physics Communications | p3[ . 

The KB-equations are designed to study time- 
dependent non-equilibrium phenomena but they can also 
be used to study the system in its final equilibrium 
state. An example is found in |^ where an initially 
un-correlated zero-temperature Fermi-distribution time- 
develops into a correlated system whilst the collisions 
are calculated with the time-evolving correlated Green's 
functions. The build-up of correlations is manifested as 
an asymptotic decrease in potential (correlation) energy 
(initially zero), with a corresponding increase in kinetic 
energy, while the total energy is constant. The resulting 
Green's functions contain a wealth of information such 
as correlation energy and particle distribution. Spectral 
functions are also easily derived. Although the collision 
term basically implies a second order calculation with 
respect to the potential the propagators are by the pro- 
cess of time-iteration dressed with second order insertions 
(with their proper energy-dependence) up to all orders. 

In this paper we further address and explore these fea- 
tures of the KB-equations. In particular we focus on 
correlation times and energies. We also extend previ- 
ous comparisons [ p^ with the Levinson equation. This 
equation can be obtained from the KB-equations by ap- 
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proximating the propagators in the coUision integral by 
free Green's functions. The comparison between the two 
approaches allow us then not only to explore the validity 
of the Levinson equation but also to asses the importance 
of the higher order diagrams associated with the dressing 
of the propagators. 

Although the interaction can in principle be a T- 
matrix, including ladder summations, the presented cal- 
culations are done with an effective time - local inter- 
action. Section 2 contains a short summary of the KB- 
formalism used in this paper. Section 3 deals with the 
Levinson equation and some relations involving correla- 
tion energies and the Born approximation are shown. An 
apparent dilemma concerning the Born approximation 
result in Section 3 is resolved in Section 4 using the ex- 
tended quasiparticle picture. Numerical results are shown 
in Section 5 and Section 6 summarizes our findings. 



II. THE KB-EQUATIONS 

We show some of the equations regarding the KB- 
formalism needed for ourpresentation. For further de- 
tails see for example refs 0J|,|l^. 

In a homogeneous medium neglecting the mean field 
the KB-equations reduce (with 1) to: 



/ dt"(S> (p, <") - E<(p, t, t"))G< (p, t", t') - 

Jto 

f dt"lf< (p, t"){G> (p, t\ t') - G< (p, t", t')) 

Jtn 



(1) 



/ dt"{G> (p, t, t") - G<(p, t, i"))S< (P, t", t') 

Jto 



dt"G< (p, t")(S> (P, i", - S<(p, 0) 



to 



(2) 



The notations are the conventional ones. and G^ 
are essentially the occupation-numbers for holes and par- 
ticles respectively. The particle distribution function 
P{P, t) is given by 



p{p,t) = -^G<(p,^,^). 



(3) 



The Green's functions G^ and G^ are related on the 
diagonal in the t, t' plane by 



G>{p,t,t) = -i + G<{p,t,t). 



(4) 



Between these two Green's functions there exists a useful 
relation 

GS(p,t,i')-[G<(p,i',t)]*. (5) 
The scattering rates S are given by 

j:<{p,t,t') = -z J |^G>(pi,i',i) 

x/^(p-Pi) |r<(p + pi,t,0 I i(p-pi)\. (6) 



Here T< is defined by 
(p\T<{P,t,t')\p 
J dt"dt"'dp"dp"' (p I T+(P,t,t") I i(p" - p'") 



xG< (p",t",r)G< (p"',i",n 
x(i(p"-p"')|T-(P,C't')|P 



(7) 



The effective interaction is usually defined in a 
binary collision (ladder) approximation by an integral 
equation formally written as 



Tf2 = V^ + VG±G±T± 



(8) 



where V is the 'free' interaction potential. 

In the following will be approximated by a local 
time-independent effective interaction 



V{p) = ^3/2^3-^^g-ir,^p^ 



(9) 



with rj = 0.57 fm and Vq = —AdSMeV. Considering the 
full dynamical T-matrix approximation one obtains more 
involved time integrals which gives rise to nonlocal effects 

The exchange term is not included. The eq. (|g) for 
the scattering rates then simplifies to 



/ /N f d'^pxd'^q „ , , 

s<(p,M') = -* J ^^^'(q) 



(27r)6 

X G> (pi , t', t)G< (q + Pi , t, t')G< (p - q, t, t'). (10) 

The momentum integrations are conveniently evaluated 
using the convolution theorem for Fourier transforms. 
The diagrammatic representation can be seen in figure 
|l|. We will compare in the paper the selfconsistent with 
the non-selfconsistent approximation which are given in 
figure by thick and thin lines respectively. 
The total energy of the system reads 



Etot{t) = ^Kp{t) + 



d^p ^ d d 



4 7 (27r)3'9ii dt2 



)G {pMM)) 
(11) 
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which is used throughout the paper for calculating the 
total energy. The kinetic energy Kp for the correlated 
medium is 



(27r)3 2m' 

and the correlation energy is defined by 
E,,„{t) = Etotit) - Kp{t). 



(12) 



(13) 



Note that the mean or Hartree-Fock field is not included 
in our work and the total energy therefore contains only 
the correlated energy. We also define the uncorrelated 
kinetic energy by 



Kf{t) 



P 



(27r)3 2m 



fip,t). 



(14) 



The relation between the reduced density matrix p and 
the quasi-particle distribution / will be discussed below. 
(See Section 4) Several numerical applications of this for- 
malism are published. For some references see the Intro- 
duction. A detailed description of numerical details are 
published in ref . [n3| . 



+ 



ns 



FIG. 1. The integral form of Kadanoff-Baym equation 
(above) and the two discussed approximations for the self- 
energy (below). The first line describes the selfconsistent 
Born approximation (thick lines) and the second line the 
non-selfconsistent Born approximation (thin lines). The lat- 
ter one leads to the Levinson equation 
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In this paper we concentrate on correlation energies 
and correlation times. All (with one exception) calcula- 
tions are performed with an initially uncorrelated nuclear 
matter system and a momentum-distribution /(p,t — 0) 
specified by a density p and temperature T. The system is 
then time-evolved beyond equilibrium. The selfenergies 
are conserving |l|Jl^ so that the total energy is conserved. 
Therefore we have during the time evolution 



(15) 



where Xj- is the kinetic (and total) energy of the initial 
unperturbed and uncorrelated system and is the to- 
tal energy after equilibration (t — > oo) including correla- 
tions. 



III. THE LEVINSON EQUATION 

As we will show below the Levinson equation can be 
considered as one of several proposed approximations of 
the KB-equations. It is of particular interest to us here 
because it allows for some analytic results [^,0 that 
help to understand the numerical findings and also help 
to illustrate some important features built into the KB- 
equations. Let us first discuss different reduction schemes 
of the KB equation including the Levinson equation. 

To evolve the two-time KB-equations along t = t' one 
needs the Green's functions for off-diagonal times t ^t' . 
Different methods have been proposed approximating the 
off-diagonal elements in terms of the diagonal ones and 
using eq.(||) a one-time theory can then be derived. 

The first method we want to discuss is the Generalized 
Kadanoff Baym (GKB)-ansatz introduced by Lipavsky 
et al. Il9] given by 



(16) 



for t' > t and 



G>(p,i,t')-G>(p,t',i')5(p,t,i') (17) 

for i > t'. S is the spectral function defined by 

5(p, t, t') = z(G> (p, <') - G< (p, i, t')). (18) 

The GKB ansatz was discussed in ref and numer- 
ical results were also shown of various approximation 
schemes. The very simplest ansatz for the spectral func- 
tion is the quasiparticle approximation to get 

G<(p,t,i') = G<(p,i,t)e-(*-*') 

G>(p,t,i') = (-Z + G<(p,t',t'))e— (19) 

The KB-cquations then reduce to Levinson's equation 
[^,^-22 1 for homogeneous systems 



^P(p,t) 



d^Pid^p'd^p'i ^ 



2 / " ''72^y "' ^(1 p - p' i)'^(p + pi - p' 



p'l) 



X / dTCOsAw(i - r)(p(p',r)p(p'i,T)p(p,T)p(pi,r) 
-P(P, t)p(pi , r)p(p', r)p(p;, r)) (20) 

with p — \ — p and Aw — (lo + lo\ — lo' ~ uj'^ . The mean 
field is neglected in the following and the quasiparticle 
energy is approximated by an effective mass w — /2m. 
If for large times the distribution function becomes sta- 
tionary the integration over the cos-function reduces to 
(5(Aa;) and the equilibrium distribution p will be a Fermi- 
distribution. 

By the approximation ([l^) the correlated (damped) 
Green's functions in the KB-equations are replaced by 
the free Green's functions which in other words means 
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that the second order energy diagram is calculated with- 
out insertions in the particle (hole) lines. (See Fig. |^). 
In the time-evolution of the KB-equations the lines are 
on the other hand dressed repeatedly with insertions. 
This results in a damping or dephasing while there is no 
damping in the Lcvinson equation but there does remain 
a memory-effect; the integration over r in eq (|2C|). 

It is instructive to notice that the Boltzmann collision 
integral is obtained from equation ( |2C| ) if: 

(i) One neglects the time retardation in the distribution 
functions, i.e. the memory effects and 

(ii) The finite initial time to is set equal to — c» corre- 
sponding to what is usually referred to as the limit of 
complete collisions. 

For the Markovian Boltzmann equation the kinetic en- 
ergy is conserved, while the potential energy is zero. The 
Levinson equation conserves the total energy. [^3| The 
correlation energy is now given by 



-^corr (^) — 

1 r d^pd^pid^p'd^p[ 



Vi\p-p' imP + Pi-p'-p'i) 



X f dTsmAuj(t -t){p{p',t)p{p[,t)p{p,t)p{pi,t) 

Jta 

-p(p, r)p(pi, t)p(p', r)p(p;, r)). (21) 
For large times eq. (^ij) reduces to ||l^ 



corr 



1 r d^pd^pid^p'd^p'^ 

2 J (2^0^ 



The second order Born approximation for the total en- 
ergy is on the other hand known from perturbation the- 
ory 



niPi-p'i I)'5(p + Pi-P'-P'i) 



Peq (Pi )Peq (P2 )Peg (Pl )Peq (P2) — 



(22) 



jpeq 
^tot 



1 



corr 



(24) 



(As noted above the Hartree-Fock energy is not included 
so that accordingly the first order contribution to the 
energy is not included here.) One should note that in 
the process of equilibration the system is excited and the 
correlation energy does in both expressions ( p3|) and ( |2^ ) 
refer to the excited but not to the initial ground state of 
nuclear matter. Also note that in the process of excita- 
tion the uncorrelated kinetic energy ify- has increased to 

f ■ 

In order to resolve the apparent disagreement between 
eqs. (23) and (HJ) we first note that eq. (|2^ ) results 
from a time-evolution of the Levinson equation starting 
from an uncorrelated system with a kinetic energy 
We shall show below that the correlated and uncorre- 
lated kinetic energies at the end of the time-evolution 
are related by 



Kf 



1 



corr 



(25) 



Upon insertion in eq. (23) the apparent disagreement is 
then resolved. 

To prove eq. (^5|) we have to discuss the difference 
between the reduced density matrix p and the quasipar- 
ticle distribution /. This is performed in the extended 
quasiparticle picture. 



IV. EXTENDED QUASIPARTICLE PICTURE 



where V denotes the principal value and where peq in- 
dicates the equilibrium large time correlated densities. 
This energy resembles the second order Born estimate of 
the potential energy but with two important differences. 

The first of these is that the densities Peq are corre- 
lated densities, in the long time equilibrated limit. A 
Born estimate would however be done with an uncorre- 
lated distribution /. For weak interactions and/or low 
density for which the Levinson equation and the Born 
approximation are certainly valid, the difference between 
initial uncorrelated and final correlated densities is neg- 
ligible. Around nuclear matter values this difference is 
however important and we shall address this question 
below showing numerical results for Levinson densities, 
KB densities as well as initially uncorrelated densities. 

The second difference between eq. (^2|) and the sec- 
ond order Born estimate is that upon closer inspection 
a factor of one-half appears missing. This will now be 
clarified. With the correlation energy given by eq. (|2^ ) 
the total energy after equilibration is using eq. (O) 



In the previous section we derived the Levinson equa- 
tion from the time diagonal parts of the Kadanoff and 
Baym equations. We adopted the GKB ansatz for the 
closure of the off-diagonal parts of the Green's functions. 
But there is another path that we now want to explore. 

Our alternative is to use the extended quasiparticle 
approximation (EQP) for the spectral function in the fi- 
nal equilibrated system. The EQP is consistent with the 
Levinson equation as both are a low density and/or weak 
interaction approximation but they differ in their phys- 
ical content and different renormalizations p4| . In an 
expansion with /mS+ << i?eE+ and dReY^^duj << 1 
one finds ||l9.Ea,E6[ for the spectral function 



Seqp{p,oj) = 2^6{uj - u:o)Z{p) - v "^^]^^^ ^"^' f (26) 



where 



K) = EZ = k;^ 



corr 



(23) 



Z(p)^l + (^^^^|^)_„ 



(27) 
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is the wave function renormalization and where V indi- 
cates that the principal value is to be taken when inte- 
grating over S. The energy is defined by 



cjo =p^/2m + ReY.+ {Tp,uJo) 



(28) 



and S+ is the retarded selfenergy. (See also Sect. 5.2). 
The EQP approximation satisfies the first two lo weighted 
sum rules 



(29) 



and has been well tested numerically for nuclear matter 

.. . . ^ . . 

In equilibrium the expansion ( |26| ) is consistent with 

the following ansatz 



G<{p,u) = z/(p)27rZ(p)(5(c^ - uo) + S<(p,c^) 



V 



{u - wo)^' 
(30) 



It allows for an approximate construction of the Green's 
function and the Wigner or reduced density matrix in 
terms of the quasiparticle distribution. The relation be- 
tween these two distribution functions, the quasiparticle 
distribution and the reduced density matrix, was first in- 
troduced by Craig |^ within the limit of small scatter- 
ing rates. An inverse functional /[p] was constructed in 
ref. [ p8| . For equilibrium nonideal plasmas the correlated 
density has been employed in refs. |p^ , p9| and under the 
name of the generalized Beth-Uhlenbeck approach it has 
been used in ref. [|30| for nuclear matter studies. The au- 
thors in refs. [ p5| , ^6[ have used this approximation under 
the name of Extended Quasiparticle Approximation for 
the study of the mean removal energy, reduced density 
and high-momenta tails of the reduced density matrix. 

The non-equilibrium (time-dependent) extension of 
this formalism has been recovered within the quasiparti- 
cle approach to the kinetic equation for weakly interact- 
ing particles and referred to as a modified Kadanoff and 
Baym ansatz |]3l| , ^. 

Integrating eq. (^0[) (in the time-dependent extension) 
over iv gives a relation between the reduced density ma- 
trix and the quasiparticle distribution 



/d 
^S<(p,L^, 



V 



(31) 



The imaginary part of the retarded function is ob- 
tained from 

/mS+(p,c^,t) ^ -^(E>(p,c^,t) - (E<(p,c^,i)) (32) 

and the real part is obtained from the dispersion relation. 
Then 



Z{Y>,t)^l-^ J J^^-^(S<(p,t)-S>(p,t)) 



27r {lo — ujq)^ 



(33) 



so that eq. (^Tj) can be rewritten as 
/ \ \ f duj V 

J ^^{UJ-UJoY 

x(E<(p,u;,t)(l - /(p,i)) + S>(p,c^,t)/(p,i)) (34) 

providing also the relation between the uncorrelated and 
correlated energies. Multiplying with the kinetic energy 
^ and integrating over p one finds the relation (|2^) that 
we wanted to prove. 



We like to point out that the first term in eq (34) is 
just the uncorrelated distribution while the second cor- 
rects for the off-shell scatterings. It was shown in ref 
psi that this equation gives practically exact agreement 
with a Brueckner calculation. One should also note that 
the Levinson equation for the reduced density-matrix was 
derived with the ansatz (|l^) i.e. with a quasi-particle 
spectral function. The correlations induced in the time- 
evolution then result in a spectral function consistent 
with the EQP-approximation. If one instead uses the 
EQP spectral function in the ansatz the result is a non- 
local Boltzmann - like kinetic equation for the quasipar- 
ticle distribution [|l6|,0. For the frequency independent 
real valued interaction used here the nonlocal effects van- 
ish and we are left with the Boltzmann equation for the 
quasiparticle distribution / from which one then can find 
the reduced density matrix by (jsj). Alternatively one 
can consider the Levinson equation for the reduced den- 
sity matrix directly. This is the approach taken in this 
paper. For more detailed discussions, derivations and 
physical content of these relations we refer the reader to 
refs 

Let us remind that the two outlined schemes, the 
Levinson equation as well as EQP, are strictly valid only 
at low density and/or for weak interactions. All the same, 
the EQP has been found to be an excellent approxima- 
tion for nuclear matter lE5| . 



V. NUMERICAL RESULTS 

In the above we have displayed the KB- as well as the 
Levinson-equations; the latter is obtained from KB us- 
ing free (quasi-particle) Green's functions in the collision 
kernel. 

In this section we show some results of our calcula- 
tions. The equations (KB and Levinson) were time- 
evolved starting at time t — Q with an uncorrelated Fermi 
gas of specified density and temperature. In addition to 
KB and Levinson correlation energies we have also cal- 
culated the second order Born energy given by eq. (|2^). 
In one case we also consider the collision between two 
Fermi-spheres (slabs in coordinate-space). 
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A comparison with the approximate Levinson and 
Born results requires a high precision of the calcula- 
tions to be meaningful. To minimize the relative er- 
rors all calculations are made with essentially the same 
KB computer-program. To perform the Levinson calcu- 
lations all that has to be changed relative to the KB is to 
replace the selfconsistently calculated Green's functions 
used in the KB-code in the collision-kernel with the free 
Green's functions shown in eq (^. In the second order 
Born-approximation calculations the Green's functions 
on the time-diagonal were replaced by constant (time- 
independent) distribution functions p(p) as specified in 
the respective calculations below. 

It may be noteworthy to point out that the Pauli- 
blocking which in most perturbative calculations such as 
Brueckner is treated approximately (e.g. with an angle 
averaging) here is treated " exactly" . 

The meshes that we use are an improvement relative to 
previous work as more computer-power is now available. 
The momentum-mesh in the Cartesian coordinate system 
that we use was 57 points along each axis with Ap^. = 
Apy = Ap^ = 0.2 fm^^. The number of time-steps was 
typically 80 with A< ~ 0.25/m/c . The interaction shown 
in eq. (^) was used. The dependence of the strength Vq 
of the interaction is now also investigated. 



A. Approach to Equilibrium; Correlation Times 

If one time-evolves an initially uncorrelated Fermi- 
distribution with the classical Boltzmann-equation (mod- 
ified for Fermion-statistics) this initial distribution would 
be stationary. The Fermi-distribution is in fact a solu- 
tion of this equation with the collision-term equal to zero. 
There are no two-body correlations. This is not so for 
the Kadanoff-Baym equation. As the system in this case 
evolves from the initial uncorrelated state the correlation 
energy decreases with time. The system correlates. At 
the same time the kinetic energy increases by the same 
amount, because our choice of self-energies is conserving 
i.e. the total energy is conserved. As a consequence of 
the self-consistent time-iterations higher and higher or- 
ders of insertion diagrams are included in the Green's 
function propagators until full selfconsistency and equili- 
bration is achieved. In the KB-case the system converges 
to an equilibrated distribution but this is not always so 
for the Levinson-case. We shall return to this below. 

In a previous publication we have already consid- 
ered the correlation time tc, i.e. the time it takes the 
system to correlate time-evolving from the initially un- 
correlated Fermi-distribution. It was in particular found 
that in the Levinson-approximation a simple result could 
be derived at low temperature namely 



tn = 



1 



(35) 



was found to be relatively small. This result was also used 
in a subsequent paper on interferometry methods |34| . 

Fig. H shows the correlation energies normalized to the 
equilibrium energy as a function of time for five different 
densities as indicated in the figure-caption. The initial 
distribution is in each case a zero temperature Fermi- 
distribution. 
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FIG. 2. Correlation energy from KB-calcuIations as a 
function of time normalized to final equilibration energies. 
From left to right the five curves correspond to the normal 
density of nuclear matter (po = O.lSfm"'^) multiplied by 
2, 1, 0.5, 0.25 and 0.125. See also Table \ 

In accordance with ref. |jl^ we define the correlation 
times tc as the time of maximum correlation-energy. We 
now find that the times tc scale roughly as 



tr 



I 



2Ef 



(36) 



with Ep the Fermi-energy. The temperature-dependence 



as shown in Table |. The previous analytic result ( |35| ) 
based on the Levinson equation gave correlation times 
twice as large. A precise definition of these times is how- 
ever not possible as is also seen in figure |^. 

Fig. 1^ shows the correlation energies normalized to the 
equilibrated energy at five different interaction-strengths 
indicated in the figure-caption. The initial distribution 
is in each case a zero-temperature Fermi-distribution. 
The previous analytic finding |lj]was that the correla- 
tion times are roughly independent of the interaction 
strengths. This is nicely confirmed by the numerical re- 
sults of Fig. H for times smaller than 2 fm/c. There is 
however a large overshoot as the strength is decreased. 
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This still remains to be understood. With the particular 
definition of the correlation-times that we have adopted 
we still see a dependence on the interaction-strengths be- 
cause of this overshoot as displayed in Table 



III 




2 4 6 8 10 12 14 
TIME(fm/c) 



FIG. 3. Similar to Fig. |^ but for normal nuclear mat- 
ter density and different interaction strengths. The normal 
strength Vb = 453 is multiplied by 2.0, f .0, 0.5.0.25 and O.f 25 
respectively. The figure shows t there is an increasing "over- 
shoot" of the correlation energy as the strength is decreased. 



See also Table III 



In a previous paper it was also shown that the temper- 
ature dependence of tc should be quite small. |Q This 
conclusion was based on the Levinson equation. The re- 
sult using the KB-equations shown in Fig. ^ confirms 
this. One may only note that the approach to equili- 
bration is markedly different at higher temperatures as 
compared to zero temperature. Our general conclusion 
regarding the independence of temperature is also ex- 
emplified by the second curve from the left in Fig. ^. 
This shows the correlation energy in a "collision" be- 
tween two Fermi-spheres separated by a momentum of 
(100 MeV/A collision energy) and with a total 
density of O.lS/m"^ (normal density). 

In the calculations with the KB-equations our choice of 
selfenergies are conserving, i.e. total energy is conserved. 
psf The Levinson-equations also conserve total energy. 
From our calculations we find after 80 time-steps a de- 
crease in total energy of 3.8% at normal nuclear matter 
density and similarly an increase of 1.3% in the KB-case. 
This numerical accuracy is quite satisfactory. 
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FIG. 4. Correlation energy as a function of time normal- 
ized to final equilibration energies. The interaction strength 
and density are both normal. From left to right the first,third 
and fourth curves are for initial temperatures 0, 20 and 40 
MeV respectively. The second curve is for a "collision" event 
(see text). Notice the near overlap of the curves. 

Another important result shown by Figs 2-4 is that 
the correlation energy (and the kinetic energy) converges 
to an equilibrium value in the KB-case. While the KB- 
results show excellent convergence this is not true in the 
Levinson case as shown by Fig. ^ which should be com- 
pared with Fig. H One finds a convergence only for the 
two weakest interactions which also agree quite well with 
the corresponding curves in Fig. |^ while only moderate 
convergence with the factor 0.5. For the strongest inter- 
action there is even sign of an instability in agreement 
with the findings of Haug |35| and others |Q , |^ who 
found a continuous increase of energy with time. The 
origin of this problem with the Levinson equation was 
discussed in terms of the artificial coUisional broadening 

We note that in the Levinson case the Green's func- 
tions are free, uncorrelated, in the collision-term while 
the ensuing Green's functions are not free as e.g. evi- 
denced by the non-zero potential energies. This is an in- 
consistency in the Levinson equation which, we believe, 
is related to the non-convergence. The problem with this 
equation is also reflected in the results in Table III show- 



ing reasonably good agreement between KB and Levin- 
son correlation energies only for the three weakest inter- 
actions. 
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FIG. 5. Time-dependence of normalized correlation ener- 
gies in the Levinson-case for normal nuclear matter density. 
The oscillating curve is with the interaction strength twice 
the normal. The smooth curves from bottom and up are with 
strengths 1., 0.5, 0.25 and 0.125 times the normal. Notice that 
the correlation energy converges only for the two weakest in- 
teractions. Compare with the KB-results in Fig. ^ The 
initial distribution is in each case a zero-temperature distri- 
bution of normal nuclear matter density. 

The fact that the correlation time is finite and of the 
order of some fm/c in nuclear matter is, we believe, of 
importance when studying heavy ion collisions. Corre- 
lations will change with density and temperature with 
typically this response time. In a previous paper we have 
already looked at the consequence of this fact on inter- 
ferometry measurements. [p4| 



B. Equilibrium Spectral Function; Occupation 
Numbers 

Correlation energies are intimately related to the spec- 
tral functions in our case given by eq. (|l8|) . In the quasi- 
particle approximation valid for weak interactions and an 
uncorrelated medium they are simply delta-functions cor- 
responding to vanishing widths represented by the first 
part of eq. (p6|). 

In this section we consider the system in its final equi- 
librated state. The total energy is then related to the 
spectral function by 



^-t-oo 



+ 00 



(27r)3^2TO 



c^]^(p,c^) (37) 



which is the same as eq. 
tions 



|11D with the equilibrium rela- 



G<(p,c.)=z5(p,o.)/H 
G>(p,c.)=z5(p,o.)(l-/H). 



(38) 



A factor 4 is included for the spin, isospin etc. degeneracy 
of the nuclear system. 

The spectral functions S'(p, w) can be calculated from 
eq. (|l^) by Fourier-transforming from the time-domain. 
It was found more practical however to use 



5(p,L^) 



2/mS+(p,w) 



2m 



i?eS+(p,a;))2 



(/mS+(p,w))2 
(39) 



Figure ^ shows some results at normal nuclear matter 
density. These are from KB-calculations. (As noted 
above in section 3 the spectral functions in the Levin- 
son approximation are in agreement with the EQP- 
approximation ( p^ ) in lowest order densities.) The ini- 
tial distribution is a zero-temperature Fermi-distribution 
and the w-dependent selfenergies are obtained by Fourier- 
transforming from the time-domain after equilibration. 
The correlation energy is in this case 35.95 MeV as shown 
m Table I The fi gure shows that the widths are compa- 
rable with those from Brueckner and other many body 
calculations with "realistic" interactions. If the calcula- 
tion were for the ground-state the distribution would be 
much more peaked at this momentum. Because the cal- 
culations in our present work are made from an initially 
uncorrelated state that is then time-evolved along the 
real axis the final state corresponds to an excited state 
(with a temperature estimated to be about 25 MeV). For 
comparison we also show in fig. ^ spectral functions for 
the ground state obtained with imaginary time-stepping. 
[^,|| The momenta at the Fermi-surface is not shown here 
as it would be too high and narrow in this case, but note 
the peaks at the adjacent momenta. 

The occupation-numbers are in equilibrium also re- 
lated to the spectral functions by 



1 /■+°° 



f{uj)S{p,uj)duj 



(40) 



consistent with eqs. (||) and (p^). 

Figure ^ shows a comparison of the density- 
distributions obtained respectively from the KB- and 
Levinson calculations. Please note that the KB calcu- 
lations lead to a pronounced discontinuity at the Fermi 
energy. A quasiparticle (Fermi-) distribution / with a 
temperature of 27 MeV is plotted for comparison. This 
is the roughly estimated temperature of the final equili- 
brated system. This would be the distribution in a Boltz- 
mann calculation in which case there are no correlations, 
the spectral function is the quasi-classical and the p and / 
distributions are identical. The KB - as well as Levinson 
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distributions possess power tails at high momenta and 
suppress states at lower momenta due to correlations [0. 

Correlation energies are shown in Tables || and 111. 
These are energies after 80 time-steps (20/to/c). At this 
time the system is well equilibrated in the case of KB as 
already seen in the figures ||, |^ and ^. In other words the 
Green's functions along the time-diagonal become con- 
stant for large times. 

This is not so for the off-diagonal elements that carry 
the correlations. Fig|| shows the absolute value of these 
elements normalized to the equilibrium value, as a func- 
tion of past times. For the KB-calculations the well- 
known damping or dephasing is seen. 





ENERGY (MeV) 



FIG. 6. Similar to Fig [7| except that these spectral func- 
tions are for the ground-state and obtained with imaginary 
time-stepping. The momenta are from left to right .0, 1.0, 1.6 
and 2.2/m"^ 

The slight increase for the largest past times is a nat- 
ural consequence of the correlation process. While the 
equilibrated occupation for p is 0.62. (See fig ^) the 
initial value was 1.0 and this is still "memorized" 20 fm/c 
later. 
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FIG. 7. Spectral functions in nuclear matter at 

p = 0.18/m~^, kf — 1.4/m~^. The momenta are from left to 
right 0.0,0.75, 1.2 and 1.65/m~^. Note that the last momen- 
tum is larger than kf. The temperature is here estimated to 
be 25 MeV. 

Of particular interest is the increase in damping as the 
strength of the interaction is increased which is well il- 
lustrated by the figure. This is a very important effect 
that is contained in the KB-equations with the selfconsis- 
tently calculated Green's functions. This has as a conse- 
quence that the memory-time decreases; the integration 
over past times which in principle should start at the 
time when interactions are switched on, {to in eqs. (|l| 
and |)), can now be limited to about 5 fm/c or less be- 
fore the point of time of the evolution. Fig 9 in ref. ||] 
illustrates this point. Associated with an increase of the 
strength of the interaction is of course not only a decrease 
in memory-time but also an increase in the width of the 
spectral functions. 

Only for the very weakest strength the damping is neg- 
ligible. This is consistent with Table |l| that shows a 
good agreement between all methods for this case. If 
the situation were such that the evolution at the same 
time is sufhciently slow another limit is reached. The 
time-integration can now be done analytically the energy- 
spread decreases and collapses to a delta-function. The 
Boltzmann limit is reached. 
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FIG. 8. Density distributions from the KB- and Levin- 
son- calculations at normal nuclear matter density. The 
KB-distribution is larger at small momenta and smaller at 
larger momenta. Notice that it shows the characteristic 
discontinuity at the Fermi-surface that is related to the 
strength- function. Compare with Fig. 1 in ref.(3) . 

The Levinson curve is smooth (it is a Fermi- 
distribution) and starts at about 0.52 at p=0. For 
comparison is also plotted an estimate of the Fermi- 
distribution / with the same nuclear matter density and 
a temperature of 27 MeV. It starts at about 0.68. 

The damping, so characteristic in the KB-case is ab- 
sent in the Levinson case. The corresponding curves can 
increase or decrease in the latter case depending on the 
past distribution- function (on the time-diagonal). In the 
present calculations this function will always decrease for 
p = and the plotted ratio in Fig. ^ is now consequently 
larger than 1.0. 



shown in columns 5,6 and 7 in table | showing Born en- 
ergies calculated with occupation numbers from the KB, 
Levinson and initial distributions (labelled "KB", "Lev." 
and "Init.") respectively. As expected the three energies 
agree exactly at the lowest density but the agreement be- 
comes progressively worse as the density is increased. As 
predicted by eq. ( p2| ) the Born-column indicated "Lev." 
(column 6) agrees on the other hand nearly exactly with 
the Levinson result (column 4) at all densities, being a 
confirmation of computational accuracy. 
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FIG. 9. This figure shows the absolute value of 
G<{p,T,t')/G<{p,T,T) for T-t' for four different strengths 
of the interaction. The four lower curves with the ratio less 
than 1.0 is from KB-calculations. From bottom to top the nor- 
mal strength is multiplied by factors of 2., 1., 0.5 and 0.125. 
The uppermost curve is from Levinson calculation with nor- 
mal interaction strength showing no damping but rather an 
increase with past times. This effect is explained in the text. 
The momentum is chosen to be p = for all cases. 



C. Equilibrium correlation energies 

It was already shown above that the Levinson corre- 
lation energy for large times approaches a second order 
Born value. But it was also pointed out with reference 
to eq. ( p^ ) that the Born calculation then has to be 
made with the Levinson final reduced density but not 
with the initial quasiparticle distribution as would nor- 
mally be done in a perturbative expansion. In the low- 
density and/or weak interaction limit, where the Born 
and Levinson approximations both become valid this dif- 
ference should be irrelevant. This is verified by the results 



1. Comparison of Levinson with Bom 

As noted above the Born(Lev) correlation energies 
agree very well with the Levinson energies at all densi- 
ties and temperatures but they differ from the Born(Init) 
energies as seen comparing columns 6 and 7 in table |. 
This is because the occupation-numbers differ in these 
two calculations. The Born(Init) calculations are made 
with the Fermi-distributions of the temperatures T in- 
dicated. These are also the initial distributions in the 
Levinson calculations, but the final distributions p (with 
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which also Born (Lev) is calculated) are different. There 
are two distinct effects contributing to the change from 
initial to final distribution. One is the heating associ- 
ated with the correlations and the second comes from 
the difference between the reduced density matrix p and 
the quasiparticle distribution / that is contained in the 
off-shell part in EQP (34). We are going to discuss these 
two effects now. 

The heating of the system in the Levinson case occurs 
since it correlates. The uncorrelated quasiparticle distri- 
bution / is consequently of a higher temperature in the 
Born(Lev) case, column 6, than in the Born(Init) case, 
column 7. 

In order to make a meaningful comparison between 
the two results we have to correct for this temperature- 
difference. Fortunately this is straightforward in the 
Levinson case. The numerical solution of the Levinson 
equation gives us the correlated distribution p{p, t) and 
from eq. ( |l2| ) the kinetic energy Kpit) in the equilibrated 
medium (for large times). The relation ( p5| ) valid for the 
Levinson equation then gives us the kinetic energy K'i^ 
which together with the density allows us to deduce the 
relevant parameters, temperature and chemical potential 
for the final equilibrium Fermi-distribution. 

The Levinson correlation energy is then compared with 
the Born(Init) correlation energy calculated with that 
same final Fermi-distribution. The result of this compar- 
ison is shown in Table || at normal density. With the 
Born and Levinson calculations compared at the same 
temperature (Tj in the table) i.e. with the same distri- 
bution f (uj) there is however a remaining difference, in 
Table || indicated by Diff. We remind that when Born 
is calculated with the p-distribution (i.e. in Born(Lev)) 
rather than with the /-distribution we get no difference, 
i.e the Levinson energy. The remaining difference Diff 
is then due to the second of the two separate effects re- 
ferred to above. It is attributed to the difference between 
the reduced density p and the corresponding quasipar- 
ticle distribution /. The former contains the spectral 
correction as in eqs. ( |3l| ) and (40). 

Such a spectral correction can be included in e.g. 
Brueckner type calculations by iteration but is rarely 
done. (See however a recent work [^). Table || shows 
a decrease of 3.3MeV of the correlation energy at zero 
temperature. This implies a decrease in binding energy 
of 1.6MeV. (We like to point out that in the EQP- 
approximation a change 6 in correlation energy for a 
given distribution f{u!) changes the total energy by 6/2. 
It changes the kinetic energy Kp by —6/2.) 

The dependence on the strength of the interaction is 
presented in Table III. The Born calculations are here 
for initial zero temperature Fermi distributions. For the 
three weakest interactions all three correlation energies 
agree reasonably well. But at full interaction strength 
we see differences of over 20 % increasing to 50 % at 
twice the strength. 



2. Comparison of KB with Levinson 

We now want to compare the KB- with the Levinson- 
equation. Although the KB collision term in our calcula- 
tions is formally up to second order in the interaction the 
correlation energy at equilibrium is actually of a much 
higher (infinite) order. This is because the correlated 
Green's functions are formed by iterative time-stepping 
functionals of the interaction. 

The effect of these higher order terms can be as- 
sessed from comparing the KB and Levinson correlation- 
energies. The difference between the two stems from the 
difference between the Green's functions in the collision 
kernel for the two separate cases. In the KB case they 
are selfconsistent (correlated) while in the latter they are 
free Green's functions. In diagrammatic language the 
presence of the correlated Green's functions in the colli- 
sion kernel means that hole- and particle- lines are dressed 
with the two-point (second order) insertions to all or- 
ders and all time orderings, see Fig. 0. This implies 
that for a similar calculation in the w-representation the 
proper dependence on uj of the selfenergy T,'^{p,uj) has 
to be included in the calculation. The same problem is 
faced if applying Brueckner 's theory in his original recipe 
where the insertions in particle lines are calculated off the 
energy-shell. This is quite complicated and 

is a compelling reason why the 'continuous' recipe (with 
both particle- and hole-insertions being on the energy- 
shell) is now customarily used. It seems that the numeri- 
cal problem with off-energy shell propagation is simplified 
by going to time-space, as in the present work. 

In Table | one finds that the difference between the 
" KB" and " Lev" correlation energies in columns 3 and 4 
at normal density amounts to 13.72 MeV for an initial 
temperature T = and decreases to 5.63 MeV at an 
initial T = QOMeV. This implies a difference in binding 
energy of 6.86 and 2.82 MeV respectively. (See discussion 
above regarding factor of 1/2.) One also finds in the same 
table that the difference decreases with density which is 
to be expected. The difference in binding energy would 
be half of the quoted numbers. The conclusion is that 
the higher order diagrams resulting from insertions in 
the propagators contained in the KB eqs but neglected 
in Levinson are quite important. 



D. Discussion 

Referring to the discussion above in table |^ our re- 
sults do actually consider excited, not ground state nu- 
clear matter. In the Levinson case we could easily es- 
timate the excitation temperature. (See section V.B). In 
the KB-case there is no such simple relation. Using the 
Levinson relation as a rough estimate we find however 
an excitation energy of about 25 MeV at normal den- 
sity. Ground state calculations can be done by imaginary 
time-stepping but this has not yet been done. The cal- 
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culation showing the zero temperature spectral functions 
in Fig. ^ is not sufficiently accurate to allow a detailed 
comparison of correlation energies. 

We do expect the relatively large difference between 
the Levinson and KB correlations to prevail at zero 
temperature. We also argue that our Levinson calcula- 
tion closely corresponds to a 'continuous' choice (on the 
energy-shell mean field insertions in hole- and particle- 
lines) Brueckner calculation. With this choice the higher 
order Brueckner diagrams that usually are considered 
(three-body, low-order ring and some fourth order) then 
contribute only a few MeV. 1 42 1 The Brueckner mean field 
insertions mentioned above are of course not present in 
the Levinson case but because we use a local interaction 
they are expected to give a negligible correction. We 
have confirmed this numerically. With realistic non-local 
(momentum-dependent) effective interactions one has of 
course important dispersion-corrections. 

With the self-consistent Green's functions as in the 
KB-calculations the width of the spectral functions is 
(implicitly) included. In contrast the Brueckner calcula- 
tions are done in a quasi-classical approximation with the 
effective interaction, the reaction-matrix, calculated as- 
suming an uncorrelated zero-temperature Fermi distribu- 
tion. Using Green's function methods normally including 
also hole-hole ladders one also rarely goes beyond the first 
quasi-classical approximation. (Except for the hole-hole 
diagrams this is equivalent then to Brueckner.) Spectral- 
functions are however readily obtained and the equations 
can be iterated. This is a rather formidable calculation, 
but has recently been done with realistic forces. The 
KB-calculations in this paper corresponds to such an iter- 
ated calculation with selconsistent spectral-functions and 
off-energy shell E's albeit with a simple interaction and 
no ladder summation. The results from the two calcu- 
lations are not readily compared but both agree in that 
the corrections going beyond the quasi-classical approxi- 
mation should be considered. 

Referring to Fig. |^ an important difference between 
the Levinson and KB-coUision terms is the damping, 
which is related to the width of the spectral-function. 
Referring to eq. ( p9| ) this is on the other hand related 
to S+(p,w). Above we deduced a 1.6 MeV difference 
in binding energy between the Born (Init) and the Born 
(Levinson) calculations stemming from a spectral correc- 
tion of the time diagonal density-matrix. The present 
difference between Levinson and KB of 6.86 MeV is also 
a spectral correction but stems from the off-diagonal el- 
ements i.e. the iv dependence of S+. 



VI. SUMMARY AND OUTLOOK 

The Kadanoff-Baym equations have been used to cal- 
culate correlations and to find the effect of approxima- 
tions such as using undressed propagators. The KB- 
calculations as presented above are fairly simple. The 



fortran version of the program is available from Com- 
puter Physics Communications fl^ . 

In our KB- as well as Levinson-calculations the ini- 
tial conditions are in general Fermi-distributions of spec- 
ified temperature and density. As the equations are time- 
stepped the system correlates with a decrease in the po- 
tential energy. The total energy is conserved and the 
kinetic energy increases. The temperature of the system 
is consequently increasing. The final correlation-energies 
are therefore not for the ground-states of the respective 
systems. To explicitly study the ground-state one can 
use the imaginary time-stepping technique. |^,^ It would 
have been desirable to use this in our present work for 
calculation of correlation energies. The precision is how- 
ever not (yet) at the level of the computations with the 
uncorrelated initial condition. Here we only used the 
imaginary time-stepping method to calculate the spec- 
tral function at zero temperature shown in Fig. |^. 

One of the purpose of this work was to investigate the 
time it takes for the system to correlate from the initially 
uncorrelated state. We have in particular calculated the 
dependence upon density, temperature and also strengths 
of the interaction as shown in Sect. 5.2. We find that 
the correlation-time tc scales roughly as 



tr. = 



2Ef 

but is nearly independent of the strength of the inter- 
action and of the temperature of the system. This time 
is relevant for the discussion of collisions between heavy 
ions when tc is comparable with the collision-time be- 
tween the ions. It is also of practical importance in 
dynamic calculations. In principle the calculation with 
the KB-equations involves an integration over all past 
times referred to as a memory-effect. In practice this is 
not necessary. The correlations effectively cuts down the 
memory-time. This is also demonstrated by the damp- 
ing shown in Fig. |^. In nuclear systems the memory-time 
is typically less than S/m/c, see ref. or 10-20 time- 
steps. In electron plasmas where correlations are smaller 
the corresponding time is 100/s or 25-50 time-steps. 

It was shown in a previous work that the Levinson 
correlations approach a second order Born expression at 
large times. |1J] Increasing the strength of the interac- 
tion an instability does however develop as shown in Fig. 
^. A dilemma involving a factor of 2 comparing with the 
Born expression for the correlation energy was resolved in 
Sect 4. In order to compare the Levinson with the Born 
correlation energies a temperature correction had to be 
applied and the result was shown in Table ^ There re- 
mains a difference between the two which amounts to a 
decrease in binding of 1.6MeV This can be ascribed to a 
spectral correction and is due the difference between the 
uncorrelated distribution /(p) and the reduced density 
p(p) in the Levinson calculation. The relation between 
/ and p is expressed by eq. (PQ) which in the EQP ap- 
proximation reduces to eq. (|34|). 
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A comparison of the KB- with the Levinson- 
correlations show the effect of the dressing of the prop- 
agators by the insertions. This is found to result in 
a comparatively large correction; a 6.85 MeV decrease 
in binding energy. This correction is also ascribed to 
a spectral correction but now for time off-diagonal ele- 
ments. This result should have an impact on the ongoing 
problem of nuclear saturation which is found to be above 
the experimental density and/or energy for perturbative 
schemes without spectral corrections. 

The total energy given by eq. includes besides a 
kinetic energy only the binary correlation energy. There 
is of course also a first order term, the Hartree-Fock term, 
contributing to the energy of the many-body system. 
This mean field term can be included easily in the Green's 
functions, jl^ We point out that besides the usual mo- 
mentum independent Hartree shift another effect appears 
usually referred to as the dispersion-effect. This effect 
stemming from the momentum-dependence of the Fock 
(or Brueckner-Hartree-Fock) field is well-known and not 
of interest in our present calculation. It further decreases 
the binding energy in nuclear many-body calculations. 

All of our nuclear matter calculations have until now 
been restricted to using a time-local interaction. The cor- 
relations appear however to be similar to those for more 
realistic interactions. This is illustrated by the spectral 
functions shown in Figs. |^ and |^. They show a width 
comparable to more serious calculations and the expected 
behavior as a function of momentum. [ p5[ Our interaction 
does not have a short-ranged repulsion or tensor-part but 
we believe the long-ranged part to be a reasonable repre- 
sentation of the true interaction. We do however envision 
a future extension to more realistic nucleon forces as well 
as a T-matrix. 

It may finally be relevant to point out that an impor- 
tant difference between the present KB-calculations and 
more conventional Green's function studies is that it has 
been performed here in t' rather than lu -space. This is 
found to be very practical. 

This work was supported in part by the National Sci- 
ence Foundation Grant No. PHY-9722050 
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TABLE L 

Correlation energies as a function of the density of nuclear matter. At normal density the temperature dependence 
is also shown. All results are here with Vq = 453.0MeV. The energies i^^orr ^'^'^ (the negative of) the equilibrium 
correlation energies. The Born energy E^^n (Born) is calculated with three different distribution-functions as discussed 
in the text. 



p 

fm 


1 

MeV 


-E-corr (KB) 

MeV 


-Bcorr (Lev) 
MeV 


KB 


EcoTT (Bo 
MeV 
Lev. 


rn) 

Init. 


tc (KB) 
fm/c 


ft 

2Ep 

fm/c 


0.380 





53.63 








65.16 


2.0 


1.5 


0.183 





35.95 


49.67 




49.69 


43.97 


2.4 


2.4 


0.181 


10 


36.03 


48.60 


50.52 


48.54 


49.16 






0.182 


20 


35.94 


46.74 




46.65 


52.84 






0.182 


40 


34.31 


42.14 




42.09 


50.52 






0.182 


60 


31.59 


37.22 


38.00 


37.20 


44.65 






0.095 





23.55 


31.33 




31.76 


28.83 


3.4 


3.8 


0.047 





14.40 


17.84 


18.24 


18.26 


17.47 


5.2 


6.2 


().()2:i 





«.12 


9.80 


l).9() 


9.9() 


l).9() 


8.5 


9.7 



TABLE IL 

Comparison of Levinson E'corr ^^'^ Born results at equal uncorrelated kinetic energies {K^^). The initial temperature 

Ti of the Levinson calculation is increased to Ty as a consequence of the correlations. Born(Init) is the Born correlation 
energy at this same temperature Tf. The remaining difference Diff. between E^^rr ^^'^ Born is due to the correlational 
spreading of the spectral function and is discussed in the text. All energies are in MeV. The density is here normal 
nuclear matter density. 



Ti 


^^c'c^.(Lev) 


Kf 




Tf 


Born(Init) 


Diff. 





49.67 


24.18 


49.03 


27 


53.0 


3.3 


10 


48.60 


29.51 


53.78 


31 


52.5 


3.9 


20 


46.74 


40.91 


64.24 


38 


51.0 


4.3 


40 


42.14 


67.93 


88.98 


54 


46.5 


4.4 


()() 


:!7.22 


9(). 10 


115.00 









TABLE IIL 



Correlation energies as a function of strength of the interaction Vq. All results are here for normal nuclear matter 
density p = 0.183fm^^ and temperature T = 0. The Born energies are here calculated only with the initial uncorrelated 
distribution. 



Vo 
MeV 


EtSrr (KB) 

MeV 


SSrr (Lev) 

MeV 


E^Srr (Born) 
MeV 


tc (KB) 
fm/c 


906 


93.15 


134.41 


175.9 


1.8 


453 


35.95 


49.67 


43.97 


2.4 


227 


10.45 


12.41 


10.99 


3.2 


113 


2.72 


2.85 


2.74 


3.5 


57 


0.68 


0.70 


0.69 


3.6 
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